Geometrical aspects of entanglement 

Jon Magne Leinaas a , Jan Myrheim c and Eirik Ovrum"' 6 
(a) Department of Physics, University of Oslo, 
P.O. Box 1048 Blindern, 0316 Oslo, Norway 
(b) Centre of Mathematics for Applications, 
P.O. Box 1053 Blindern, 0316 Oslo, Norway and 
(c) Department of Physics, The Norwegian University 
of Science and Technology, 7491 Trondheim, Norway 
(Dated: May 8, 2006) 

Abstract 

We study geometrical aspects of entanglement, with the Hilbert-Schmidt norm defining the metric on 
the set of density matrices. We focus first on the simplest case of two two-level systems and show that a 
"relativistic" formulation leads to a complete analysis of the question of separability. Our approach is based 
on Schmidt decomposition of density matrices for a composite system and non-unitary transformations to a 
standard form. The positivity of the density matrices is crucial for the method to work. A similar approach 
works to some extent in higher dimensions, but is a less powerful tool. We further present a numerical 
method for examining separability, and illustrate the method by a numerical study of bound entanglement 
in a composite system of two three-level systems. 

PACS numbers: 03.67. -a, 03.65.Ud, 02.40.Ft 
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INTRODUCTION 



Entanglement is considered to be one of the main signatures of quantum mechanics, and in 
recent years the study of different aspects of entanglement has received much attention. One ap- 
proach has been to study the formal, geometrical characterization of entangled states as opposed 
to non-entangled or separable states [1-3]. In such a geometrical approach the Hilbert-Schmidt 
metric defines in many ways the natural metric on the space of of physical states. This metric 
follows quite naturally from the Hilbert space norm, when the quantum description is extended 
from pure to mixed states, and it is a Euclidean metric on the set of density matrices. For a com- 
posite quantum system the separable states form a convex subset of the full convex set of density 
matrices, and one of the aims of the geometrical approach is to give a complete specification of 
this set and thereby of the non-separable or entangled states. 

The purpose of this paper is to examine some questions related to the geometrical description 
of entanglement. We focus primarily on the simplest composite systems consisting of two two- 
level systems (2x2 system) or two three-level systems (3x3 system), but examine also some 
questions relevant for higher dimensions. In the case of two two-level systems the separable states 
can be fully identified by use of the Peres criterion [4]. This criterion states that every separable 
density matrix is mapped into a positive semidefinite matrix by partial transposition, i.e., by a 
transposition relative to one of the subsystems. Since also hermiticity and trace normalization is 
preserved under this operation, the partial transpose of a separable density matrix is a new density 
matrix. 

A non-separable density matrix, on the other hand, may or may not satisfy Peres' condition. 
This means that, in general, Peres' condition is necessary but not sufficient for separability. How- 
ever, for the special case of a 2 x 2 system as well as for a 2 x 3 system the connection is stronger, 
and the Peres condition is both necessary and sufficient for a density matrix to be separable [5]. 
To identify the separable density matrices, and thereby the entangled states, is therefore in these 
cases relatively simple. In higher dimensions, in particular for the 3 x 3 system, that is not the 
case. Peres' condition is there not sufficient for separability, and the convex subset consisting of 
all the density matrices that satisfy this condition is larger than the set of separable matrices. This 
is part of the reason that the identification of the set of separable states in higher dimensions is a 
hard problem [6]. 

In the present paper we first give, in Sect. 2, a brief introduction to the geometry of density 
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matrices and separable states. As a next step, in Sect. 3 we focus on the geometry of the two- 
level system and discuss the natural extension of the Euclidean three dimensional Hilbert-Schmidt 
metric to a four dimensional indefinite Lorentz metric. This indefinite metric is useful for the 
discussion of entanglement in the 2x2 system, where Lorentz transformations can be used to 
transform any density matrix to a diagonal standard form in a way which preserves separability 
and the Peres condition (Sect. 4). By using this standard form it is straightforward to demonstrate 
the known fact that any matrix that satisfies the Peres condition is also separable. 

The transformation to the diagonal standard form is based on an extension of the Schmidt de- 
composition to the matrix space with indefinite metric. For general matrices such a decomposition 
cannot be done, but for density matrices it is possible due to the positivity condition. We show 
in Sect. 5 that the Schmidt decomposition in this form can be performed not only for the 2x2 
system, but for bipartite systems of arbitrary dimensions. However, only for the 2 x 2 system 
the decomposition can be used to bring the matrices to a diagonal form where separability can be 
easily demonstrated. In higher dimensions the Schmidt decomposition gives only a partial simpli- 
fication. This indicates that to study separability in higher dimensions one eventually has to rely 
on the use of numerical methods [7-1 1]. 

In Sect. 6 we discuss a new numerical method to determine separability [12], and use the 
method to study entanglement in the 3 x 3 system. The method is based on a numerical estimation 
of the Hilbert-Schmidt distance from any chosen density matrix to the closest separable one. We 
focus particularly on states that are non-separable but still satisfy the Peres condition. These are 
states that one usually associates with bound entanglement. A particular example of such states is 
the one-parameter set discussed by P. Horodecki [13] and extended by Bruss and Peres [14]. We 
study numerically the density matrices in the neighbourhood of one particular Horodecki state and 
provide a map of the states in a two dimensional subspace. 

THE GEOMETRY OF DENSITY MATRICES 
The convex set of density matrices 

A density matrix p of a quantum system has the following properties, 

pt = p hermiticity 
p > positivity 
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Tr p = 1 normalization (1) 
The matrices that satisfy these conditions form a convex set, and can be written in the form 

P = ^2pk\^k)(^k\ (2) 

k 

with \ip k ) as Hilbert space unit vectors and 

Pk>0, 5> fe = l (3) 

k 

The coefficient p k , for given k, can be interpreted as the probability that the quantum system is in 
the pure state \^ k ). This interpretation depends, however, on the representation (2), which is by 
no means unique. In particular, the vectors \ip k ) may be chosen to be orthonormal, then they are 
eigenvectors of p with eigenvalues p k , and Eq. (2) is called the spectral representation of p. 

The pure states, represented by the one dimensional projections are the extremal points 

of the convex set of density matrices. That is, they generate all other density matrices, correspond- 
ing to mixed states, by convex combinations of the form (2), but cannot themselves be expressed 
as non-trivial convex combinations of other density matrices. 

The hermitian matrices form a real vector space with a natural scalar product, Tr(AB), which is 
bilinear in the two matrices A and B and is positive definite. From this scalar product a Euclidean 
metric is derived which is called the Hilbert-Schmidt (or Frobenius) metric on the matrix space, 

ds 2 = \ Tr((dp) 2 ) (4) 

The scalar product between pure state density matrices pi = \ip){tp\ and p 2 = \4>){4>\ is 

Tr( PlP2 ) = |<V#>| 2 (5) 

For an infinitesimal displacement \dtp) on the unit sphere in Hilbert space the displacement in the 
matrix space is 

dp=|#)M + M<#| (6) 

and the Hilbert-Schmidt metric is 

ds 2 = (#|#) - |M#)| 2 (7) 

where we have used that Tr(rfp) = {ip\dip) + {dip\ip) = 0. This may be interpreted as a metric 
on the complex projective space, called the Fubini-Study metric. It is derived from the Hilbert 
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space metric and is a natural measure of distance between pure quantum states, in fact, ds is the 
infinitesimal angle between rays (one dimensional subspaces) in Hilbert space. Since the Hilbert- 
Schmidt metric on density matrices is a direct extension of the Hilbert space metric, it is a natural 
metric for all, both pure and mixed, quantum states. 

A complete set of basis vectors {J a }, in the space of hermitian matrices, can be chosen to 
satisfy the normalization condition 



For n x n matrices the dimension of the matrix space, and the number of basis vectors, is n 2 . One 
basis vector, J , can be chosen to be proportional to the identity 1, then the other basis vectors are 
traceless matrices. A general density matrix can be expanded in the given basis as 



where the coefficients £ a are real, and the trace normalization of p fixes the value of £ - With the 
chosen normalization of J we have £ = l/v^-- 

Due to the normalization, the density matrices are restricted to a hyperplane of dimension 
n 2 — 1, shifted in the direction of J relative to a plane through the origin. The set of density 
matrices is further restricted by the positivity condition, so it forms a closed, convex set centered 
around the point p = J /y/n. This point corresponds to the maximally mixed state, which has the 
same probability (ip\p \ip) = 1/n for any pure state The geometry is schematically shown in 
Fig. 1, where the set of density matrices is pictured as the interior of a circle. One should note that 
the normalization condition in a sense is trivial and can always be corrected for by a simple scale 
factor. In the discussion to follow we will find it sometimes convenient to give up this constraint. 
The quantum states can then be viewed as rays in the matrix space, and the positivity condition 
restricts these to a convex sector (the cone in Fig. 1). 

Unitary transformations 

The hermitian matrices J a can be viewed as generators of unitary transformations, 



Tr (J a Jb) = Sab 



(8) 




(9) 



a 




(10) 
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FIG. 1 : A schematic representation of the set of density matrices in the vector space of hermitian matrices. 
The positive matrices form a cone about the axis defined by the unit matrix, and the normalization condition 
restricts the density matrices to a convex subset, here represented by the shaded circle. 

with real coefficients ( a , which act on the density matrices in the following way, 

p ^~p = XjpXj\ (11) 

If p is represented as in Eq. (2), then we see that 

p = ^2pk\i>k)(i>k\ (12) 
k 

where \-ipk) = U\ipk)- Thus, the matrix transformation p — > UpW is induced by the vector 
transformation — > U\ip). An immediate consequence of Eq. (12) is that the transformed 
density matrix UpW is positive. 

Such unitary transformations respect both the trace and positivity conditions and therefore leave 
the set of density matrices invariant. Also the von Neumann entropy 

S = -Tr(p\np) (13) 

is unchanged, which means that the degree of mixing, as measured by the entropy, does not change. 

Since the dimension of the set of density matrices grows quadratically with the Hilbert space 
dimension n the geometry rapidly gets difficult to visualize as n increases. However, the high 
degree of symmetry under unitary transformations simplifies the picture. The unitary transforma- 
tions define an SU (n) subgroup of the rotations in the n 2 — 1 dimensional matrix space, and all 
density matrices can be obtained from the diagonal ones by these transformations. In this sense 
the geometry of the set of density matrices is determined by the geometry of the set of diagonal 
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density matrices. The diagonal matrices form a convex set with a maximal set of n commuting 
pure states as extremal points. Geometrically, this set is a regular hyperpyramid, a simplex, of 
dimension n — 1 with the pure states as corners. The geometrical object corresponding to the full 
set of density matrices is generated from this by the SU (n) transformations. 

In Fig. 2 the set of diagonal density matrices is illustrated for n = 2, 3 and 4, where in the 
first case the hyperpyramid has collapsed to a line segment, for n = 3 it is an equilateral triangle 
and for n = 4 it is a tetrahedron. For n — 2, the SU (2) transformations generate from the line 
segment the three dimensional Block sphere of density matrices. This case is special in the sense 
that the pure states form the complete surface of the set of density matrices. This does not happen 
in higher dimensions. In fact, the dimension of the set of pure states is 2n — 2, the dimension of 
SU(n)/(U(l) x SU(n — 1)), because one given pure state has a U(l) x SU(n — 1) invariance 
group. This dimension grows linearly with n, while the dimension of the surface, n 2 — 2, grows 
quadratically. 




n^2 n=3 n=4 



FIG. 2: Geometric representation of the diagonal density matrices for the three cases of Hilbert space 
dimension 1, 2 and 3. For general dimension n they define a hyperpyramid of dimension n — 1. 

The faces of the hyperpyramid of dimension n — 1 are hyperpyramids of dimension n — 2, cor- 
responding to density matrices of the subspace orthogonal to the pure state of the opposite corner. 
Similarly, the hyperpyramid of dimension n — 2 is bounded by hyperpyramids of dimension n — 3, 
etc. This hierarchy is present also in the full set of density matrices, generated from the diagonal 
ones by SU(n) transformations. Thus, to each extremal point (pure state) the boundary surface 
opposite to it is a flat face corresponding to the set of density matrices of one lower dimension. In 
this way the boundary surface of the set of density matrices contains a hierarchy of sets of density 
matrices of lower dimensions. 
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The boundary of the set of density matrices is characterized by at least one of the eigenvalues 
of the density matrices being zero, since outside the boundary the positivity condition is broken. 
This means that at the boundary the density matrices satisfy the condition det p = 0, which is 
an algebraic equation for the coordinates of the boundary points. When n is not too large the 
equation can be solved numerically. This has been done in Fig. 6 where a two dimensional section 
of the set of density matrices is shown. One should note that there will be solutions to the equation 
det p = also outside the set of density matrices. The boundary of the set of density matrices can 
be identified as the closed surface, defined by det p = 0, that encloses the maximally mixed state 
and is closest to this point. 

More general transformations 

We shall later make use of the complex extension of the transformations (10), by allowing ( a to 
be complex. This means that the transformation group is extended from SU(n) to SL(n, C) (the 
normalization condition det U = 1, or ( = 0, is trivial). Transformations of the form p = VpV^ 
do not respect the trace condition Tr p = 1 if V is non-unitary, but they do respect the positivity 
condition, because they are still vector transformations of the form (12), with \ip k ) = V\^k)- This 
means that they leave the sector of non-normalized density matrices invariant. They no longer 
keep the entropy unchanged, however. Thus, the larger group SL(n, <D) connects a larger set of 
density matrices than the restricted group SU (n) . 

One further generalization is possible. In fact, even if we allow V to be antilinear, the transfor- 
mation p = VpV^ still preserves positivity, because Eq. (12) still holds with \ip k ) = V\ip k ). This 
point needs some elaboration. 

An operator V is antilinear if 

V(a\i>) + b\(f))) = a*V\4>) + b*V\<P) (14) 

for any vectors \4>) and complex numbers a, b. Let {\i)} be a set of orthonormal basis vectors, 
let ipi = {i\ip) and write \Vip) for the vector V Then 

\ViP) = Vj2^\i)=J2 V r^\ti W 
with Vji = (j\V\i). The hermitian conjugate is defined in a basis independent way by the 
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identity 

Mvty) = (16) 

or equivalently, 

|yV) = yt|0) = ^^..| i ) (17) 

By definition, V is antiunitary if V* = V -1 . 

Familiar relations valid for linear operators are not always valid for antilinear operators. For 
example, when V is antilinear and \ Vip) = V\ip), it is no longer true that {Vip\ = {ip\V^. This 
relation cannot hold, simply because (Vip\ is a linear functional on the Hilbert space, whereas 
{ip\V^ is an antilinear functional. What is nevertheless true is that V\^){t/j\V^ = \Vip)(Vip\ . In 
fact, both of these operators are linear, and they act on the vector \<f>) = £\ 4>j\j) as follows, 

v\^){^\4>) = v [y^^v^m ) = ^iV&vw) = \v^)(vm m 

\ i,j J i,j 

As a consequence of this identity the form (12) is valid for the antiunitary transformations, and the 
positivity is thus preserved. 

The transposition of matrices, p — > p T , obviously preserves positivity, since it preserves the set 
of eigenvalues. This is not an SU(n) transformation of the form (11), as one can easily check. 
However, transposition of a hermitian matrix is the same as complex conjugation of the matrix, 
and if we introduce the complex conjugation operator K, which is antilinear and antiunitary, we 
may write 

p T = KpK^ (19) 

Note that transposition is a basis dependent operation. The complex conjugation operator K is also 
basis dependent, it is defined to be antilinear and to leave the basis vectors invariant, K \i) = 
We see that = K = K~ x . 

One may ask the general question, which are the transformations that preserve positivity of 
hermitian matrices. If we consider an invertible linear transformation on the real vector space 
of matrices, then it has to be a one-to-one mapping of the extremal points of the convex set of 
positive matrices onto the extremal points. In other words, it is a one-to-one mapping of one 
dimensional projections onto one dimensional projections. In yet other words, it is an invertible 
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vector transformation \ip) — > V\ip), defined up to a phase factor, or more generally an arbitrary 
non-zero complex factor, for each pure state One can show that these complex factors can be 
chosen in such a way that V becomes either linear or antilinear, and that the matrix transformation 
is p — > VpVK However, we will not go into details about this point here. 

Geometry and separability 

We consider next a composite system with two subsystems A and B, of dimensions ua and ub 
respectively. By definition, the separable states of the system are described by density matrices 
that can be written in the form 

P = ^Pkp£®pj! (20) 
k 

where p^ and pf are density matrices of the two subsystems and pk is a probability distribution 
over the set of product density matrices labelled by k. The separable states form a convex subset of 
the set of all density matrices of the composite system, with the pure product states where 
\ip) = |0) <8> |x), as extremal points. Our interest is to study the geometry of this set, and thereby 
the geometry of the set of entangled states, defined as the complement of the set of separable states 
within the full set of density matrices. 

The Peres criterion [4] gives a necessary condition for a density matrix to be separable. Let us 
introduce orthonormal basis vectors \i)a in Ha and \ j) B in Hb, as well as the product vectors 

lu) = V)a ® \j)b (21) 
We write the matrix elements of the density matrix p as 

Pij-M = (ij\p\kl) (22) 
The partial transposition with respect to the B system is defined as the transformation 

P^ P P ■ Pij-M = Pil;kj (23) 

This operation preserves the trace, but not necessarily the positivity of p. However, for separable 
states one can see from the expansion (20) that it preserves positivity, because it is just a transpo- 
sition of the density matrices pf of the subsystem B. 
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Thus, the Peres criterion states that preservation of positivity under a partial transposition is a 
necessary condition for a density matrix p to be separable. Conversely, if the partial transpose p p 
is not positive definite, it follows that p is non-separable or entangled. The opposite is not true: if 
p p is positive, the density matrix p is not necessarily separable. 

It should be emphasized that the Peres condition, i.e., positivity of both p and p p , is indepen- 
dent of the choice of basis vectors \i) A and \j)b- In fact, a change of basis may result in another 
definition of the partial transpose p p , which differs from the first one by a unitary transforma- 
tion, but this does not change the eigenvalue spectrum of p p . The condition is also the same if 
transposition is defined with respect to subsystem A rather than B. This is obvious, since partial 
transposition with respect to the A subsystem is just the combined transformation p — > (p T ) p ■ 

Let us consider the Peres condition from a geometrical point of view. We first consider the 
transposition of matrices, p — > p T , and note that it leaves the Hilbert-Schmidt metric invariant. 
Being its own inverse, transposition is an inversion in the space of density matrices, or a rotation 
if the number of inverted directions, n(n — l)/2 in an n x n hermitian matrix, is even. Since p and 
p T have the same set of eigenvalues, transposition preserves positivity and maps the set of density 
matrices onto itself. Thus, the set of density matrices V is invariant under transposition as well as 
under unitary transformations. 

Similarly, a partial transposition p — > p p preserves the metric and therefore also corresponds 
to an inversion or rotation in the space of matrices. On the other hand, it does not preserve the 
eigenvalues and therefore in general does not preserve positivity. This means that the set of density 
matrices, V, is not invariant under partial transposition, but is mapped into an inverted or rotated 
copy V p . These two sets will partly overlap, in particular, they will overlap in a neighbourhood 
around the maximally mixed state, since this particular state is invariant under partial transposition. 
We note that, even though partial transposition is basis dependent, the set of transposed matrices 
V p does not depend on the chosen basis. Nor does it depend on whether partial transposition is 
defined with respect to subsystem A or B. 

To sum up the situation, we consider the following three convex sets. V is the full set of density 
matrices of the composite system, while V — V n V p is the subset of density matrices that satisfy 
the Peres condition, and S is the set of separable density matrices. In general we then have the 
following inclusions, 

S C V C V (24) 



11 



The Peres criterion is useful thanks to the remarkable fact that partial transposition does not 
preserve positivity. This fact is indeed remarkable, for the following reason. We have seen that 
any linear or antilinear vector transformation \ip) — > V\ip) will preserve the positivity of hermitian 
matrices by the transformation p — > VpVK It would seem that V — 1a <S> K b , a complex 
conjugation on subsystem B, would be a vector transformation such that p p = V pV\ and hence, 
that partial transposition would preserve positivity. What is wrong with this argument is that there 
exists no such operator as 1a <E> K b . To see why, choose a complex number c with |c| = 1, and 
consider the transformation of a product vector = \<j>) ® \x), 

(1 A ® K B M = (1 A ® K B )((c |0)) ® (c* = c 2 (|0) ® (K\ X ))) (25) 

The arbitrary phase factor c 2 invalidates the attempted definition. 

The boundary of the congruent (or reflected) image V p of V is determined by the condition 
det p p = in the same way as the boundary of the set of density matrices V is determined by 
det p = 0. As a consequence, to determine whether a density matrix p belongs to the set V is not a 
hard problem. One simply checks whether the determinants of p and p p are both positive for every 
p on the line segment between p and the maximally mixed state p . However, to check whether a 
density matrix is separable and thus belongs to the subset S is in general not easy, even though the 
definition (20) of separability has a simple form. The exceptional cases are the systems of Hilbert 
space dimensions 2 x 2, 2 x 3 or 3 x 2, where S = V. 

Schmidt decomposition and transformation to a standard form 

A general density matrix of the composite system can be expanded as 

a=0 b=0 

where the coefficients are real, and and J p are orthonormal basis vectors of the two sub- 
systems. We may use our convention that Jq = 1a/ ' \Jua~ and J p = 1 B / \/n~B~, and that and 
J p with k > are generators of SU(ua) and SU (n B ). 

A Schmidt decomposition is a diagonalization of the above expansion. By a suitable choice of 
basis vectors and J p , depending on p, we may always write 

n\-l 
a=0 
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assuming that n A < n B . There exist many different such diagonal representations of a given p, in 
fact it is possible to impose various extra conditions on the new basis vectors. It is usual to impose 
an orthonormality condition, that the new basis vectors should be orthonormal with respect to 
some positive definite scalar product. Then the Schmidt decomposition of p is the same as the 
singular value decomposition [15] of the n 2 A x n B matrix £ ab . Below, we will introduce a Schmidt 
decomposistion based on other types of extra conditions. 

The usefulness of the representation (27) is limited by the fact that we expand in basis vectors 
depending on p. However, we may make a transformation of the form 

p p = VpV ] (28) 

where V — Va <8> V B is composed of transformations Va £ SL(ua, <D) and V B £ SL(n B , <D) that 
act independently on the two subsystems and transform the basis vectors and into VaJ^a 
and V B J^V B . A transformation of this form obviously preserves the set S of separable states, 
since a sum of the form (20) is transformed into a sum of the same form. It also preserves the set 
V of density matrices satisfying the Peres condition. In fact, it preserves the positivity not only of 
p, but also of the partial transpose p p , since 

(p) P = (V A ® V* b )p p {Va ® V^ (29) 

Here V B is the complex conjugate of V B . What is not preserved by the transformation is the trace, 
but this can easily be corrected by introducing a normalization factor. Such transformations have 
been considered e.g. by Cen et al. [16]. 

As we will later show, it is possible to choose the transformation V — Va <S> V b in such a way 
that the transformed and normalized density matrix p can be brought into the special form 

/ n\-l 



1 



n A n B 



V 



k=l 



with and as new sets of traceless orthonormal basis vectors. We have some freedom in 
choosing Va and V B , because the form (30) is preserved by unitary transformations of the form 

U = U A ® U B . 

In the case n A = n B = 2 we may choose V A and V B so that j£ and are fixed sets of basis 
vectors, independent of the density matrix p. In particular, we may use the standard Pauli matrices 
as basis vectors. In this way we define a special form of the density matrices p, which we refer 
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to as the standard form. Any density matrix can be brought into this form by a transformation 
that preserves separability and the Peres condition. All matrices of the resulting standard form 
commute and can be simultaneously diagonalized. This makes it easy to prove the equality S = V, 
and thereby solve the separability problem. Although this result is well known, the proof given 
here is simpler than the original proof. 

The decomposition (30) is generally valid, but when either n A , n B , or both are larger than 2, it 
is impossible to choose both j£ and to be independent of p. Simply by counting the number 
of parameters one easily demonstrates that this cannot be done in higher dimensions. Thus, the 
product transformations V A <8> Vb are specified by 2n 2 A + 2n 2 B — 4 parameters, while the number of 
parameters in (30) is n\ — 1, when the generators j£ and are fixed. This gives a total number 
of parameters 3n 2 A + 2n 2 B — 5, when n A < n B , compared to the number of parameters of the 
general density matrix, which is n 2 A n 2 B — 1. Only for n A = n B = 2 do these numbers match. 

The mismatch in the number of parameters shows that the independent transformations V A 
and Vb performed on the two subsystems are less efficient in simplifying the form of the density 
matrices in higher dimensions. In particular, it is impossible to transform all density matrices to 
a standard form of commuting matrices. Thus, the question of separability is no longer trivially 
solved. Nevertheless, we consider the Schmidt decomposition to be interesting and important, if 
only because the number of dimensions in the problem is reduced. We expect this to be useful, 
even if, at the end, separability can only be determined by numerical methods. 

THE TWO-LEVEL SYSTEM 

The density matrices of a two-level system describe the states of a qubit, and represent a simple, 
but important, special case. It is well known that the normalized density matrices, expressed in 
terms of the Pauli matrices <x = (<7i, <7 2 , 0-3) as 

p=i(l + r-<r) (31) 

can geometrically be pictured as the interior of a three dimensional unit sphere, the Block sphere, 
with each point identified by a vector r. The two-level case is special in that the pure states form 
the complete surface |r| = 1 of the set of density matrices. The diagonal 2x2 matrices, in any 
chosen basis, correspond to a line segment through the origin with the two pure basis states as end 
points. 
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The two-level system is special also in the sense that the Euclidean metric of the three dimen- 
sional space of density matrices can be extended in a natural way to an indefinite metric in four 
dimensions. The extension is analogous to the extension from the three dimensional Euclidean 
metric to the four dimensional Lorentz metric in special relativity. Since it is useful for the discus- 
sion of entanglement in the two qubit system, we shall briefly discussed it here. 

We write the density matrix in relativistic notation as 

P=\^a lt (32) 

where a is the identity matrix 1. The trace normalization condition is that Tr p = x° = 1. We may 
relax this condition and retain only the positivity condition on p, which means that a; is positive 
and dominates the vector part r of the four-vector x M , as expressed by the covariant condition 

4 detp = (a; ) 2 - |r| 2 = x% = g^afx" > (33) 

In other words, the four-vector x^ is restricted by the positivity condition to be either a time-like 
vector inside the forward light cone, or a light-like vector on the forward light cone. The light-like 
vectors correspond to pure states, and the time-like vectors to mixed states. As already discussed, 
all points on a given line through the origin represent the same normalized density matrix (see 
Fig. 3). 







p > 
\\Qet p = { 

Tr p =1 




— K X TTJ 


>— 







FIG. 3: A relativistic view of the density matrices of a two-level system. The positive hermitian matrices 
are represented by the forward light cone and a density matrix by a time-like or light-like ray (red line). 
The standard normalization Tr p = x° = 1 breaks relativistic invariance (green horizontal line), but may be 
replaced by the relativistically invariant normalization det p = = 1/4 (green hyperbola). 
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Positivity is conserved by matrix transformations of the form 

p ^~ p = VpV \ = \ U ^ Gv (34) 

If we restrict V to be linear (not antilinear), invertible, and normalized by det V = 1, then it 
belongs to the group SL(2, <D), and L is a continuous Lorentz transformation (continuous in the 
sense that it can be obtained as a product of small transformations). Thus, preservation of positivity 
by the SL(2, <D) transformations corresponds to preservation of the forward light cone by the 
continuous Lorentz transformations. 

In order to compare the indefinite Lorentz metric to the Euclidean scalar product Tr(AB) in- 
troduced earlier, we introduce the operation of space inversion, 

a li = {l,a)^ = -a) (35) 

It is obtained as a combination of matrix transposition, or equivalently complex conjugation, which 
for the standard Pauli matrices inverts the sign of a y = a 2 , and a rotation of n about the y-axis. 
Thus, for a general hermitian 2x2 matrix A it acts as 

A -> A = RA T B) (36) 

where 

/o-i\ 

R=\ (37) 

We may also write A T = KAK\ where K = is complex conjugation. 
The Lorentz metric is now expressed by 

Trfaa,,) = 2 9lu , (38) 

and the Lorentz invariant scalar product between two hermitian matrices A = a^a^ and B = b^a^ 
is 

Tr(AB) = 2a"b fi (39) 

The invariance of the scalar product (38) can be seen directly from the transformation properties 
under SL{2, <D) transformations, 

A -> A' = VAV ] A^ A' = V^AV' 1 (40) 
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Thus, A and A transform under contragredient representations of SX(2, <D). 
The Lorentz transformed Pauli matrices 

^ = Vojfi = L\a v , V E SL(2, C) (41) 

satisfy the same metric condition (38) as the standard Pauli matrices a^. Conversely, any set 
of matrices satisfying relations of the form (38) are related to the Pauli matrices by a Lorentz 
transformation, which need not, however, be restricted to the continuous part of the Lorentz group, 
but may include space inversion or time reversal, or both. 

It is clear from the relativistic representation that any density matrix can be reached from the 
maximally mixed state p = 1/2 by a transformation corresponding to a boost. The Lorentz boosts 
generate from p a three dimensional hyperbolic surface (a "mass shell") where x^x^ = 1. This 
surface will intersect any time-like line once and only once. Thus, any mixed state is obtained 
from p by a unique boost. However, the pure states, corresponding to light-like vectors, can only 
be reached asymptotically, when the boost velocity approaches the speed of light, here set equal 
to 1. The form of the SL(2, C) transformation corresponding to a pure boost is 

f(*)=expQ€-ff) (42) 

where the three dimensional real vector £ is the boost parameter, called rapidity. Since the boost 
matrices are hermitian, a density matrix defined by a boost of the maximally mixed state will have 
the form 

P = N(S)V(S) 2 = N(£) (cosh f 1 + sinh f £ • a) , £ = | (43) 

where iV(£) is a normalization factor determined by the trace normalization condition. The nor- 
malized density matrix is 

p= i(l + tanh££-<r) (44) 

Thus, the boost parameter £ gives a representation which is an alternative to the Bloch sphere 
representation. The relation between the parameters r and £ is that 

r = tanh£| (45) 

which means that r can be identified as the velocity of the boost, in the relativistic picture, with 
|r| = 1 corresponding to the speed of light. We note that the positivity condition gives no restric- 
tion on £, and the extremal points, i.e., the pure states, are points at infinity in the £ variable. 
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ENTANGLEMENT IN THE 2x2 SYSTEM 



We consider now in some detail entanglement between two two-level systems. We will show 
that with the use of non-unitary transformations the density matrices can be written in a standard- 
ized Schmidt decomposed form. In this form the question of separability is easily determined and 
the equality of the two sets V and S is readily demonstrated. 



Schmidt decomposition by Lorentz transformations 

We consider transformations V — Va® Vb composed of SL(2, <D) transformations acting inde- 
pendently on the two subsystems, and therefore respecting the product form (20) of the separable 
matrices. We will show that by means of such transformations any density matrix of the composite 
system can be transformed to the form 

p = \ + dk ° k ® ° k j (46) 

which we refer to as the standard form. Note that the real coefficients dk must be allowed to take 
both positive and negative values. 

We start by writing a general density matrix in the form 

p = cT <g> a v (47) 

The transformation V = Va <E> V b produces independent Lorentz transformation L A and L B on 
the two subsystems, 

p _, p = V P V ] = c^Op ® cj v (48) 



with 



= L%L Bc c^ (49) 



The Schmidt decomposition consists in choosing La and L B in such a way that c^ v becomes 
diagonal. We will show that this is always possible when p is strictly positive. 

Note that the standard Schmidt decomposition, also called the singular value decomposition, 
involves a compact group of rotations, leaving invariant a positive definite scalar product. The 
present case is different, because it involves a non-compact group of Lorentz-transformations, 
leaving invariant an indefinite scalar product. 
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The positivity condition on p plays an essential part in the proof. It states that 

M/#> > (50) 

where \ip) is an arbitrary state vector. Let us consider a density matrix p which is strictly positive 
so that (50) is satisfied with > and not only with >, and let us restrict \ip) to be of product form, 
= \<f>) <8> \x) • With p expressed by (47), the positivity condition then implies 

c^m^n u > (51) 

with 

"fy = » n v = (xWAx) (52) 

These two four- vectors are on the forward light cone, in fact, it is easy to show that 

m M m M = n M n M = , m = n = 1 (53) 

We note that by varying the state vectors |0) and \x) all directions on the light cone can be reached. 
The inequality (51) holds for forward time-like vectors as well, because any such vector may be 
written as a linear combination of two forward light-like vectors, with positive coefficients. We 
may actually write a stronger inequality 

c^m^ >C>0 (54) 

valid for all time-like or light- like vectors m and n with m — n — 1. In fact, this inequality 
holds because the set of such pairs of four-vectors (m, n) is compact. 
Now define the function 

f(m,n) = — (55) 
y/m p m p n a n a 

It is constant for m and n lying on two fixed, one dimensional rays inside the forward light cone. 
It goes to infinity when either morn becomes light-like, because 

ft \ >s Cm n 

f(m, n)>— (56) 

Using again a compactness argument, we conclude that there exist four-vectors m and n such that 
f(m, n) is minimal. We may now choose the Lorentz transformations La and L B such that 

L % = ™m » L V = (57) 
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assuming the normalization conditions m^m^ = n^n p = 1. This defines La and L B uniquely up 
to arbitrary three dimensional rotations. Define 

f(m,n) = — E= (58) 
y/mfrrip n a n a 

Since f(m,n) = f(m,h), with fn p = L P A m p and n M = L p Bfi n p , it follows that f(m,n) has a 
minimum at m = n = (1, 0, 0, 0). The condition for an extremum at m — n — (1, 0, 0, 0) is that 
c ok = c k0 = for k = 1,2, 3, so that 

3 3 
k=l 1=1 

The coefficient c 00 is the minimum value of f(m, n), and hence positive. 

The last term of Eq. (59) can be diagonalized by a standard Schmidt decomposition, and by a 
further normalization p can be brought into the form (30). Finally, a unitary transformation p — > 
UpW of the product form U = Ua® Ub may be performed, where the unitary matrices Ua and 
U b may be chosen so that U A J^U\ = an d U bJ^b = This is aways possible, because 
£77(2) transformations generate the full three dimensional rotation group, excluding inversions. 
In this way we obtain the standard form (46). 

Note that the standard form (46) of a given density matrix p is not unique, because there exists 
a discrete subgroup of 24 unitary transformations that transform one matrix of this form into other 
matrices of the same form. This group includes all permutations of the three basis vectors a k ® a k , 
as well as simultaneous reversals of any two of the basis vectors. It is the full symmetry group of 
a regular tetrahedron. If we want to make the standard form unique we may, for example, impose 
the conditions d\ > d 2 > \d 3 \, allowing both positive and negative values of d 3 . 



Density matrices in the standard form 

The density matrices of the standard form (46) define a convex subset of lower dimension than 
the full set of density matrices. It is a three dimensional section of the 15 dimensional set of density 
matrices, consisting of commuting (simultaneously diagonalizable) matrices. The eigenvalues, as 
functions of the parameters d k of Eq. (46), are 

Ai, 2 = \ (1 ± (di - d 2 ) + d 3 ) , A 3 , 4 = \ (1 ± (di + d 2 ) - d 3 ) (60) 
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The pure states (with eigenvalues 1, 0, 0, 0) that are the extremal points of the convex set of com- 
muting matrices are specified by the conditions 

\di\ = \di\ = \ds\ = 1 , didzds = — 1 (61) 

There are four such states, corresponding to the four corners of the tetrahedron of diagonal density 
matrices, and these are readily identified as four orthogonal Bell states (maximally entangled pure 
states). 

We now consider the action of the partial transposition on the tetrahedron of diagonal den- 
sity matrices. It preserves the standard form, and transforms the coefficients as d\ — > d\ , d 2 — > 
— c? 2 , dz — > d$. Thus it produces a mirror image of the tetrahedron, by a reflection in the d 2 direc- 
tion (see Fig. 4). The density matrices of standard form belonging to the set V, i.e., satisfying the 
Peres condition that they remain positive after the partial transposition, form an octahedron which 
is the intersection of the two tetrahedra. 




FIG. 4: Geometric representation of the diagonal matrices spanned by orthogonal Bell states of the 2x2 
system. The density matrices form a tetrahedron (green) while the matrices obtained from this by a partial 
transposition defines a mirror image (blue). The separable states are defined by the intersection between the 
two tetrahedra, and forms an octahedron (in the center). 

We will now show that for the density matrices of standard form the Peres condition is both 
necessary and sufficient for separability. What we have to show is that all the density matrices 
of the octahedron are separable. Since, in general, the separable matrices form a convex set, it is 
sufficient to show that the corners of the octahedron correspond to separable states. 
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The density matrices of the octahedron satisfy a single inequality 

Kl + \d 2 \ + \d 3 \ < 1 (62) 

and its six corners are (di, d 2 , d 3 ) = (±1, 0, 0), (0, ±1,0), (0, 0, ±1), corresponding to the mid- 
points of the six edges of each of the two tetrahedra. The corners are separable by the identities 

- (1 ± a k <g> a k ) = ^ [(1 + a k ) <g> (1 ± a k ) + (1 - a k ) a k )\ (63) 

This completes our proof that the Peres condition is both necessary and sufficient for separabil- 
ity of density matrices on the standard form. Furthermore, since any (non- singular) density matrix 
can be obtained from a density matrix of standard form by a transformation that preserves both 
separability and the Peres condition, this reproduces the known result that for the 2 x 2 system 
the set of density matrices that remain positive after a partial transposition is identical to the set of 
separable density matrices. 

With this we conclude the discussion of the two-qubit system. The main point has been to 
show the usefulness of applying the non-unitary Lorentz transformations in the discussion of sep- 
arability. Also in higher dimensions such transformations can be applied in the form of SL(n, C) 
transformations, although not in precisely the same form as with two-level systems. 



HIGHER DIMENSIONS 



The relativistic formulation is specific for the two-level system, but some elements can be 
generalized to higher dimensions. We consider first a single system with Hilbert space dimension 
n. Again, if the trace condition is relaxed, the symmetry group SU (n) of the set of density matrices 
is extended to SL(n, C). The Hilbert-Schmidt metric is not invariant under this larger group, but 
the determinant is invariant for any n. However, it is only for n = 2 that the determinant is 
quadratic and can be interpreted as defining an invariant indefinite metric. 

The generalization to n > 2 of the Lorentz boosts are the hermitian matrices 

V = , det V = 1 (64) 
and expressed in terms of the SU (n) generators J they have the form 

V = exp(£ • J) (65) 
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with real group parameters Here J and £ are n 2 — 1 dimensional vectors. Any strictly positive 
density matrix p (with det p > 0) can be factorized in terms of hermitian matrices (65) as 

p = NV 2 = N exp(2£ • J) (66) 

with the normalization factor 

" = ttexpU • J) = (67) 
Thus, in the same way as for n — 2, any strictly positive density matrix can be generated from 
the maximally mixed state by an SL(n, C) transformation of the form (65). The boundary matri- 
ces, however, which satisfy det p = 0, cannot be expressed in this way, they can be reached by 
hermitian transformations only asymptotically, as |£| — > oo. In this limit N — > 0, and therefore 
Tr V 2 — > oo for the non-normalized density matrix V 2 . 



Schmidt decomposition 

We now consider a composite system, consisting of two subsystems of dimension ua and n B , 
and assume p to be a strictly positive density matrix on the Hilbert space H = Ha <8> T~Lb of 
the composite system. The general expansion of p in terms of the SU (n) generators is given by 
(26). Our objective is to show that by a transformation of the form (28), with V A e SL(ua, C) 
and Vg £ SL^b,^), followed by normalization, the density matrix can be transformed to the 
simpler form (30). 

Let T>a and V B be the sets of density matrices of the two subsystems A and B. The Cartesian 
product T>a xT> B , consisting of all product density matrices pa® Pb with normalization Tr pa = 
Tr p B = 1, is a compact set of matrices on the full Hilbert space H. For the given density matrix 
p we define the following function of p A and p B , which does not depend on the normalizations of 

p A and p B , 

f(n n \ Tr(p ( Pa ®Pb)) 

/(PA ' PB) ~ (det PA )Vn. (det PB y/n B 

This function is well defined on the interior of T>a x V b , where det pa > and det p B > 0. 
Because p is assumed to be strictly positive, we have the strict inequality 

Tt(p( P a® p B )) > (69) 
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and since V A x V b is compact, we have an even stronger inequality on V A xV B , 

Tr(p(p A ®p B ))>C>0 (70) 

with a lower bound C depending on p. It follows that / — > oo on the boundary of xV B , where 
either det p A = or det p# = 0. It follows further that / has a positive minimum on the interior 
of T>a xPg, with the minimum value attained for at least one product density matrix r A <g> r B with 
det t a > and det r B > 0. For t a and we may use the representation (66), written as 

t a = T\T a , t b = T b T b (71) 

ignoring normalization factors. The matrices T A e SL(n A: <D) and T B G SL(n B ,<C) may be 
chosen to be hermitian, but they need not be, since they may be multiplied from the left by arbitrary 
unitary matrices. We further write T — T A <8> T B , so that 

r A ®r B = T f T (72) 

Now define a transformed density matrix 

p = TpT* (73) 

and define 



= (detp^V- (det PB )Vn s = f(T A PAT A , T B p B T B ) 



(74) 



This transformed function / has a minimum for 

Pa ® Ps = (T^- 1 ^ <g> r^T- 1 = 1 A (g) 1 B = 1 (75) 
Since / is stationary under infinitesimal variations about the minimum, it follows that 

Ti{p5{p A ®p B )) = (76) 

for all infinitesimal variations 

8(p A ® p B ) = 8p A ®t B + l A ®5p B (77) 
subject to the constraints det (1,4 + 8p A ) = det (1b + 5p B ) = 1, or equivalently, 

Tr(5p A ) = Tr(5p B ) = (78) 
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The variations satisfying the constraints are the general linear combinations of the SU generators, 

5 PA = Jf , $PB = Y1 6c f J f ( 79 ) 



It follows that 



Tr(p Jf ® 1 B ) = Tr(pl A ® Jf) = (80) 



for all SU (n A ) generators Jf and all SU (ns) generators Jf. This means that the terms propor- 
tional to Jf <g> t B and 1 A <S> Jf vanish in the expansion for p, which therefore has the form 



fc=l i=l 



In order to write p in the Schmidt decomposed form (30), we have to make a change of basis, from 
the fixed basis sets jf and Jf to other orthonormal SU generators jf and Jf depending on p. 
This final Schmidt decomposition involves a standard singular value decomposition of the matrix 
^ki by orthogonal transformations. We may make further unitary transformations U = U A ® U B , 
but as already pointed out, this is in general not sufficient to obtain a standard form independent 
of p. 



NUMERICAL APPROACH TO THE STUDY OF SEPARABILITY 

In higher dimensions the Peres condition is not sufficient to identify the separable states. In 
other words, there exist entangled states that remain positive after a partial transposition. This is 
known not only from general theoretical considerations [5], but also from explicit examples [13]. 
States of this type have been referred to as having bound entanglement. However, whereas it is a 
fairly simple task to check the Peres condition, it is in general difficult to identify the separable 
states [6]. 

In this section we discuss a general numerical method for identifying separability, previously 
introduced in [12]. It is based on an iterative scheme for calculating the closest separable state and 
the distance to it, given an arbitrary density matrix (test state). The method can be used to test 
single density matrices for separability or to make a systematic search to identify the boundary of 
the set of separable density matrices. After giving an outline of the method we show how to apply 
the method in a numerical study of bound entanglement in a 3 x 3 system. 
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Outline of the method 



Assume a test state p has been chosen. This may typically be close to the boundary of the set of 
states that satisfy the Peres condition. Let p s be a separable state, an approximation in the iterative 
scheme to the closest separable state. We may start for example with p s = t/n, the maximally 
mixed state, or with any pure product state. The direction from p s to p is denoted a = p — p s . 
In order to improve the estimate p s we look for a pure product state p p that maximizes the scalar 
product 



or equivalently, maximizes s' = Tr(p p cr) (see Fig. 5). If s > 0, then it is possible to find a closer 



FIG. 5: Schematic illustration of the method of finding the separable state closest to a test matrix p. The 
matrix p s represents the best approximation to the closest separable state at a given step in the iteration, 
while p.p is the product matrix that maximizes the scalar product with the matrix a = p — p s . This matrix 
is used to improve p s at the next step in the iteration. The colored area S in the figure illustrates a section 
through the set of separable density matrices. 

separable state p' s by mixing in the product state p v . This search for closer separable states is 
iterated, either until no pure product state p p can be found such that s > 0, which means that 
p s is already the unique separable state closest to p, or until some other convergence criterion is 
satisfied. 

There are two separate mathematical subproblems that have to be solved numerically in this 
scheme. The first problem is to find the pure product state maximizing the scalar product s'. The 
second problem is the so called quadratic programming problem: given a finite number of pure 
product states, to find the convex combination of these which is closest to the given state p. Our 



s = Tr((p p - p a )(T) 



(82) 



S 
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approach to these two problems is described briefly below. We refer to reference [12] for more 
details. 

To approach the first subproblem, note that a pure product state p p has matrix elements of the 
form 

(ij\p p \kl) = faXj<f>* k Xi (83) 
where J2i I0«| 2 = 12 j \Xj\ 2 — 1- We want to find complex coefficients fa and Xj that maximize 

s' = Tr( Pp a) = J2 <i>i X*j ^im fa Xi (84) 

The following iteration scheme turns out in practice to be an efficient numerical method. It may 
not necessarily give a global maximum, but at least it gives a useful local maximum that may 
depend on a randomly chosen starting point. 

The method is based on the observation that the maximum value of s' is actually the maximal 
eigenvalue p in the two linked eigenvalue problems 

^2 A ik <Pk = pfa , B i l Xl = ^ ( 85 ) 

k i 

where 

A ik = ^2 x *j ai r> kl Xl ' B H = S $ ai r> kl ^ k ^ 

3,1 

Thus, we may start with any arbitrary unit vector \x) = HjXj \j)b £ an d compute the 
hermitian matrix A. We compute the unit vector |0) = fa \i) A e Ha as an eigenvector of A 
with maximal eigenvalue, and we use it to compute the hermitian matrix B. Next, we compute 
a new unit vector \x) as an eigenvector of B with maximal eigenvalue, and we iterate the whole 
procedure. 

This iteration scheme is guaranteed to produce a non-decreasing sequence of function values 
s', which must converge to a maximum value p,. This is at least a local maximum, and there 
corresponds to it at least one product vector \ fa) <g> \x) and product density matrix p p = (|0) (fa) ® 

(IxXxD- 

The above construction of p p implies, if s > 0, that there exist separable states 

p' s = (l- X)p s + Xp P (87) 
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with < A < 1, closer to p than p s is. However, it turns out to be very inefficient to search only 
along the line segment from p s to p p for a better approximation to p. It is much more efficient to 
append the new p p to a list of product states p pk found in previous iterations, and then minimize 



which is a quadratic function of coefficients \ k > with J2k — 1- We solve this quadratic 
programming problem by an adaptation of the conjugate gradient method, and we throw away a 
given product matrix p pk if and only if the corresponding coefficient becomes zero when F 
is minimized. In practice, this means that we may construct altogether several hundred or even 
several thousand product states, but only a limited number of those, typically less than 100 in the 
cases we have studied, are actually included in the final approximation p s . 

Bound entanglement in the 3x3 system 

For the 3x3 system (composed of two three-level systems) there are explicit examples of 
entangled states that remain positive under a partial transposition. This was first discussed by 
Horodecki [13] and then an extended set of states was found by Bruss and Peres [14]. We apply 
the method outlined above to density matrices limited to a two dimensional planar section of the 
full set. The section is chosen to contain one of the Horodecki states and a Bell state in addition to 
the maximally mixed state, and the method is used to identify the boundary of the separable states 
in this two dimensional plane. 

Since the separable states S are contained in the set V of states that remain positive under 
partial transposition, we start the search for the boundary of S with a series of states located at 
the boundary of V. This boundary is found by solving the algebraic equations det p = and 
det p p = 0. For each chosen state we find the distance to the closest separable state and change 
the test state p on a straight line between this point and the maximally mixed state, in a step of 
length equal to the evaluated distance. In a small number of steps the intersection of the straight 
line and the boundary of the separable states is found within a small error, typically chosen to be 
10~ 6 . (The distance from the maximally mixed state to the pure states in this case is d = 2/3.) 

In Fig. 6 we show a plot of the results of the calculations. The numerically determined points 
on the border of the set of separable states S are indicated by black dots, while the border of the set 
of states V that satisfy Peres' condition, determined by solving the algebraic equations, is shown 




(88) 
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FIG. 6: The boundary of the set of separable states in a two dimensional section, as determined by the 
numerical method. S denotes the set of separable states, B the states with bound entanglement and £ the 
entangled states violating the Peres condition. The blue straight line is determined by the algebraic equation 
det p = and gives the boundary of the full set of density matrices. The red curves are determined by the 
equation det p p = 0, in particular, the red straight line gives the boundary of the set of states that satisfy 
the Peres condition. The black dots are the numerically determined points on the boundary of the set of 
separable states. Note that the coordinates given in this plot differs from the distance defined in the main 
text by a factor y/2. 

as blue and red lines which cross at the position of the Horodecki state. One should note that the 
states with bound entanglement in the corner of the set V cover a rather small area. 

CONCLUSIONS 

To summarize, we have in this paper focussed on some basic questions concerning the geometry 
of separability. The simplest case of 2 x 2 matrices has been used to demonstrate the usefulness 
of relaxing the normalization requirement Tr p = 1. Thus, if this condition is replaced by det p = 
1, a relativistic description with a Minkowski metric can be used, where all (non-pure) states 
can be connected by Lorentz transformations. For a composite system consisting of two two- 
level systems, independent Lorentz transformations performed on the two sub-systems can be 
used to diagonalize an arbitrary density matrix in a way that respects separability. We have used 
this diagonalization to demonstrate the known fact that the Peres condition and the separability 
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condition are equivalent in this case. 

Although the diagonalization with Lorentz transformations is restricted to the composite 2 x 
2 system, we have shown that the generalized form of the Schmidt decomposition used in this 
diagonalization can be extended to higher dimensions. The decomposition involves the use of 
non-unitary SL(n, C) transformations for the two subsystems. Although a full diagonalization is 
not obtained in this way, we suggest that the Schmidt decomposed form may be of interest in the 
study of separability and bound entanglement. 

A third part of the paper has been focussed on the use of a numerical method to study sepa- 
rability. This method exploits the fact that the set of separable states is convex, and is based on 
an iterative scheme to find the closest separable state for an arbitrary density matrices. We have 
demonstrated the use of this method in a numerical study of bound entanglement in the case of a 
3x3 system. A further study of separability with this method is under way. 
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